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The effects of microwave radiation on transport properties of atomically thin La^-xSrxCuOi films 
were studied in the 0.1-13 GHz frequency range. Resistance changes induced by microwaves were 
investigated at different temperatures near the superconducting transition. The nonlinear response 
decreases by several orders of magnitude within a few GHz of a cutoff frequency u cu t ~ 2 GHz. 
Numerical simulations that assume ac response to follow the dc V-I characteristics of the films re¬ 
produce well the low frequency behavior, but fail above v cu t . The results indicate that 2D supercon¬ 
ductivity is resilient to high-frequency microwave radiation, because vortex-antivortex dissociation 
is dramatically suppressed in 2D superconducting condensate oscillating at high frequencies. 


Transport properties of thin superconducting films 
have attracted much interest due to a fascinating physical 
phenomenon, the Berezinskii-Kosterlitz-Thouless (BKT) 
transition, predicted to occur at the critical temperature 
T c lower than the temperature of the superconducting 
transition in bulk samples QJd]. The BKT phase tran¬ 
sition originates from long-range (logarithmic) interac¬ 
tions between vortex excitations in a two-dimensional 
(2D) superconducting condensate. Below T c , the dom¬ 
inant thermal excitations are vortex-antivortex (V-AV) 
pairs. A superconducting current does not move pairs, 
thus producing no energy dissipation. However, the cur¬ 
rent can break some of V-AV pairs, generate free vor¬ 
tices, set them in motion via the Lorentz force, and thus 
make the transport dissipative m\- The induced disso¬ 
ciation of V-AV pairs depends on the current strength 
and thus results in an extraordinary violation of the 
Ohm’s law. With this motivation, the strong nonlinear 
response to the electric current in thin superconducting 
films was investigated extensively [EHm* Despite signif¬ 
icant progress, important many-body and edge effects in 
the V-AV pair dissociation are still under debate mm- 
Note that the majority of studies of the nonlinear trans¬ 
port in the BKT regime have in fact been done in the 
dc domain. How condensate oscillations affect the V- 
AV dissociation is still unclear. Finally, due to strong 
phase fluctuations the BKT phenomena are significantly 
enhanced in superconducting cuprates, especially in het¬ 
erostructures with a few superconducting copper oxide 
layers EH- Our recent studies of dc nonlinearities in 
MBE-grown heterostructures show that the vortex non¬ 
linearity in the low resistive state exceeds the heating 
nonlinearity by up to four orders in magnitude m- 

Here, we present experimental investigations of 
nonlinear transport properties of atomically thin 


La 2 - x Sr x CuOi films, over a broad frequency range from 
dc to 13 GHz. The experiments indicate a dramatic de¬ 
crease of the nonlinear response at high drive frequencies, 
suggesting significant reduction of the V-AV pair disso¬ 
ciation in the oscillating superconducting condensate. 

The experiments were performed on thin 
La 2 - x Sr x Cu 04 : (LSCO) films synthesized by Atornic- 
Layer-by-Layer Molecular Beam Epitaxy, providing 
precise atomically thin layers [SHU]- On the extreme 
level of control, delta-doping in a single CuO 2 layer has 
been demonstrated P3- Recently a linear ac response 
of such films in the BKT state has been studied OS mi- 
present samples have three distinct layers. The top 
and the bottom layers, each 5 unit cells (UC) thick, 
are made of strongly overdoped (x =0.41) normal 
metals. The sample A, with 5 UC thick inner layer of 
Lai.nSrQ ^sCuOi, shows the BKT transition at T c « 
7K [21]. The sample B with 1.5 UC thick inner layer 
of Lai. 8 oSro. 2 oCuC >4 has T c « 5K. Below we study 
nonlinear response of the BKT state at T > T c . 

The films were patterned into the shape of Hall-bar 
devices with the width W =200 p,m and the distance be¬ 
tween the voltage contacts L = 800 /iin. A direct current 
I was applied through a pair of current contacts, and 
the longitudinal dc voltage V was measured between the 
potential contacts. The sample and a calibrated ther¬ 
mometer were mounted on a cold copper finger in vac¬ 
uum. The electromagnetic (EM) radiation was guided 
by a rigid coaxial line and applied to samples as shown 
in the upper insert in Fig. |T] A 50 resistor terminates 
the end of the coax and provides the broadband match¬ 
ing of the EM circuit. The EM power P and amplitude 
of the microwave voltage V v at the end of the coax were 
measured in-situ using the non-selective bolometric re¬ 
sponse of the 50 f l resistor [21]. In what follows, the 
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FIG. 1: Dependence of the voltage V on current I at different 
EM power applied to 50 D resistor at frequency v = 370 MHz. 
The thick solid line represents the dc response with no radi¬ 
ation applied. The open circles are the results of numerical 
simulations of the effect of EM radiation on I 7-1 curve com¬ 
puted at the same power, for the contact resistance R c = 210 
± 30 12 in series with the sample B at T = 7 K. The upper 
inset: a schematic of dc circuit isolated from EM circuitry by 
a capacitor C = 47 pF and three low-pass filters (LPF, v < 
100 MHz). The lower inset shows sheet resistance per square 
as a function of temperature. 


measured nonlinear response is normalized with respect 
to the calibrated EM power P, which takes into account 
all the effects of EM transmission (reflection) to (from) 
the sample stage. 

Fig. 0 shows the dependence of the voltage V on the 
current I taken at different powers P of a low frequency 
radiation. The thick solid line presents the V-I depen¬ 
dence with no radiation applied. The EM radiation in¬ 
creases the resistance as shown by thin solid lines. To 
evaluate numerically the radiation effect the electrical 
connection (coupling) between the coax and the sample 
was approximated by a high-frequency contact resistance 
R c , which determines the total voltage Vs(t) and current 
Is{t) applied to the sample, assuming that the electro¬ 
magnetic response follows the dc nonlinear V-I depen¬ 
dence [2T]. The time averaged < Vs{t) > and < Is(t) > 
are denoted by the open circles in Fig. [l] The resistance 
R c was used as the single fitting parameter for each com¬ 
puted curve, providing a good agreement between the 
low frequency experiments and the simulations. 

At small currents, both the experiment and the sim¬ 
ulations display a linear relation between the current I 
and voltage V - Ohm’s Law. Fig. [2 Ip presents the de¬ 
pendence of the ohmic resistance on microwave (MW) 
power (v = 1.5 GHz) taken at different temperatures. 
The filled dots are the slopes of the V-I dependencies at 
small currents. The power dependence varies consider¬ 
ably with the temperature. Close to the superconduct- 



FIG. 2: (a) Dependence of resistance on radiation power 
shows strikingly different behavior at low frequency (370 
MHz) and high frequencies. At frequency 3 GHz the power 
has been linearly scaled down by factor 2 to emphasize a func¬ 
tional difference in the response to EM radiation. The thick 
solid line is the numerical simulation at contact resistance 
R c ~ 210 12. Dashed line presents expected dependence at a 
small power, (b) Power dependence of resistance at different 
temperatures as labeled, v = 1.5 GHz. Sample B. 


ing state the microwave-induced resistance variations are 
strong, whereas near the normal state the variations are 
weak. The shape of the dependence changes with the 
temperature. In particular, at T = 6 K the dependence 
looks more linear than at T = 6.5 K. 

Significant changes of the power dependence are found 
in the response to microwaves with different frequencies. 
Fig. [§i shows the power dependencies of the resistance 
obtained at frequencies v = 0.37, 1.5, and 3 GHz. At low 
frequency v — 0.37 GHz, the dependence is in a good 
agreement with the numerical simulation. At high fre¬ 
quencies, the nonlinear response is much weaker and has 
a different functional form. To highlight the difference, 
the values of MW power at frequency 3 GHz were scaled 
down by a factor of 2. A comparison of power dependen¬ 
cies indicates that at a low power the nonlinear response 
is considerably weaker at a high frequency (3 GHz) than 
at a low frequency (0.37 GHz). At a high power the 
strength of the high frequency nonlinearity is restored. 

Fig. [3]shows the dependence of radiation-induced vari¬ 
ations of the dc voltage AF on dc bias I taken at different 
frequencies and power levels, as labeled. At small cur¬ 
rents, the response is linear, indicating no observable ac 
rectification in the device. The left panel demonstrates 
the effect of the low-frequency radiation (v = 0.575 GHz) 
on the resistance. All three curves are in very good agree¬ 
ment with the ones obtained by numerical simulations for 
the same radiation powers [2T| . The simulations replicate 
all the details of the experiments including the shift of 
the observed maximum with increased power using single 
fitting parameter R c = 63 ± 2 H. The right panel shows 
the effect of the high frequency (v = 3.4 GHz) radiation. 
One can see that the high-frequency response does not 
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FIG. 3: (Color online). The change of voltage AV = V(P) — 
1/(0) due to applied EM radiation as a function of the dc 
bias I at different microwave powers as labeled. The solid 
lines are experimental data; the open circles present numerical 
simulations at R c = 63 ± 2Q. Comparison of the left and 
right panels shows the agreement between experiment and 
simulation at low frequencies, but serious disparity at high 
frequencies. Sample A. 


follow the dc V-1 curve and cannot be explained by a 
reduction (re-scaling) of the MW current. 

At small currents, the linear part of the voltage vari¬ 
ations AV (/) ~ / is related to the change of the ohmic 
resistance: A R = R(P) — R( 0) = A1/(/)//. Figure[4]dis- 
plays the behavior of the resistance variation A R with 
the power at different frequencies. At a small power, 
the induced resistance variations are proportional to the 
power. At higher powers, the dependence is weaker. At 
higher frequencies, the transition to weak power depen¬ 
dence occurs at a higher power. At the highest power 
levels all dependencies converge (see Fig. [2^). 

An analysis indicates several distinct features of the dc 
nonlinear response presented in Fig. [T] At small currents 
the V - I dependence is well approximated by a combi¬ 
nation of linear and cubic terms m ■ The dependence is 
presented below: 

V = Rq ■ I + y/ 3 , (1) 

where i?o is the Ohmic resistance, the coefficient 7 is 
a constant. The high-current behavior is in agreement 
with the one expected within the BKT scenario: V ~ I a 
EM- The exponent a(T) decreases from 8 to 1 as 
the temperature increases, indicating a BKT transition 
at a(T c ) = 3 f2T]. In accordance with Eq.Q at small ac 
currents I v (t), the voltage variations A V(I) =< V(I + 
I v (t)) > —V(I) ss 3 • /• < /„(f ) 2 > are proportional to 
the dc bias /, which agrees with Fig. [3j and proportional 
to the RF power P ~ I 2 , which agrees with Fig. El 
The decrease of the nonlinear response shown in Fig. [3] 



FIG. 4: (Color online) The power dependence of radiation 
induced resistance variations A R = R(P) — R( 0) = AVjI 
obtained at small dc biases AV ~ / (see Fig. |3}. The solid 
line is the dependence A R ~ P expected for small power P. 
Open circles are numerical simulations at R c = 62 Q. Sample 
A. 


at high dc biases and the transition to a weaker power 
dependence presented in Fig. [4]are the results of the high 
current response V ~ I a at a < 2 . 

Fig. [5]presents the frequency dependence of the nonlin¬ 
ear response, taken at different temperatures and small 
powers P. At low temperatures, 9 K < T <11 K, the 
response decreases by 3-4 orders in magnitude, as the 
frequency increases above few GHz. This is in a good 
agreement with the reduction of the dc nonlinearity by 
3-4 order of magnitude observed between the regimes of 
V-AV depairing at 10 K and electron heating in the nor¬ 
mal state (see Fig. 3b in |l5] l. 

To be sure that the observed effect is not related to 
a strong decrease in the microwave coupling, we have 
evaluated the MW current through samples by investi¬ 
gating the MW reflection in the same setup (5Tj. These 
studies indicate some frequency dispersion in the MW 
current. However, the dispersion is nearly uniform and 
significantly smaller (a variation by a factor of 3-4) than 
the observed reduction of the nonlinear response by 3-4 
orders of magnitude. Independent measurements of MW 
voltage, current, and the nonlinear response in several 
samples allow us to conclude that the observed reduc¬ 
tion of the nonlinearity is of fundamental nature. 

Since the reduction is strong, we associate it with 
the frequency suppression of the dominant mechanism 
in BKT regime - the current induced V-AV pair dissoci¬ 
ation. If in the course of high-frequency oscillation the 
distance between the vortex and antivortex within a pair 
does not exceed a critical distance l c , then the pair can 
survive. At higher power, the amplitude of V-AV oscil¬ 
lations may exceed the critical distance l c , making the 
response similar to the one obtained at low frequency. It 
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FIG. 5: Frequency dependence of the relative variations of the 
resistance A R/R normalized by applied power P at different 
temperatures as labeled. Dashed line is an approximation of 
the nonlinear response between 3-6 GHz. Sample A. 


corresponds to Fig. S* and Fig. g} 

An analysis of the oscillating vortex motion inside a 
V-AV pair in the presence of dc depairing current shows 
rich physics, analogous to the Kapitza pendulum with a 
vibrating pivot point [221 . Position x v of a vortex with 
the mass M moving in a medium with the viscosity r) un¬ 
der effect of a time-dependent external potential U (r, t) 
is given by 


M 


d 2 x v 
dt 2 


dx v 


dU(r , t) 

dr ’ 


( 2 ) 


The potential U is a sum of V-AV interaction and a 
potential of Lorentz force 0 : U(r,t) = 2 qlln{r/0 - 
qoJ s (t)r, where £ is a size of the vortex core and q 2 = 
7m s 7i 2 /4TO, where n s and 2m are density and the mass 
of superconducting carriers. We express the supercur¬ 
rent density J s = Jdc + Jui as a sum of the dc bias and 
current oscillations at angular frequency to = 2ttv. 
An analysis of ac vortex motion for critical pairs near 
the saddle point of the potential ( dU/dr ss 0) shows 
a reduction of the vortex displacement at high frequen¬ 
cies: x“ c ~ qoJu/t—Mcu 2 + iurj). Moreover, at w > 
w c = 77 /M, the ac depairing is strongly suppressed due 
to the phase shift between the vortex displacement and 
the barrier height. At the moments when the potential 
U(t) reaches a minimum, the displacement is the small¬ 
est, thus preventing the V-AV dissociation. Finally at 
high frequencies a nonlinear analysis [22J of Eq. |2| indi¬ 
cates an effective attractive force inside the pair, which is 
proportional to the square of the vortex displacements: 
5F eff ~ < (z“ c (f )) 2 >~ <(Xl r ; i ) > . These effects 

reduce the pair dissociation at high driving frequencies. 
Although the origin of the cutoff frequency v cu t requires 
further research, we note that an evaluation of the fre¬ 


quency v c = is in a good agreement with our 

data. Using the Stephen-Bardeen viscosity [23] and con¬ 
sidering the vortex mass M as a total mass of electrons 
in the core [ 23 ] we get the characteristic frequency of 
i/ c = 2 GHz at the core radius of 8 nrn, which is ap¬ 
proximately two times bigger than the superconducting 
coherence length £. The above analysis points toward 
the crucial importance of the displacement-force phase 
relations in the nonlinear ac response of the bulk BKT 
state supporting substantially the fundamental origin of 
the observed phenomenon. 

In summary, strong nonlinear response to low fre¬ 
quency radiation is observed in atomically thin super¬ 
conducting films, in BKT state. The response decreases 
by several decades within a few GHz above the cutoff fre¬ 
quency v cut ss 2 GHz. This indicates that 2D supercon¬ 
ductivity is quite resilent to the high frequency radiation 
because of a strong reduction of the vortex-antivortex 
dissociation in oscillating 2D superconducting systems. 
This general conclusion is in agreement with the results 
of the linear response studies of the BKT state in thin 
disordered films of traditional superconductors. In par¬ 
ticular, detailed investigations of the conductivity in the 
range 9-120 GHz show only frequency-dependent Drude 
absorption without any measurable dissipation due to 
vortex-antivortex dissociation [25]. 

Sample synthesis by ALL-MBE and characterization 
(I.B.) and device fabrication (A.T.B.) were supported 
by U.S. Department of Energy, Basic Energy Sciences, 
Materials Sciences and Engineering Division. Work at 
CCNY and SUNY was supported by National Science 
Foundation, Division of Electrical, Communications and 
Cyber Systems (ECCS 1128459). 

SUPPLEMENTAL MATERIAL. 

TRANSPORT IN DC DOMAIN. BKT 
TRANSITION. 

According to the BKT theory, V-I characteristics 
should display a V ~ I a dependence where the 
temperature-dependent exponent, a(T), ranges from one 
in the normal state and increases as temperure decreases, 
passing through it’s critical value a(T c ) = 3 at the 
Brezenskii-Kosterlitz-Thouless (BKT) transition tem¬ 
perature, T c im ia ns]. Figure presents V-I depen¬ 
dencies of Sample A plotted in log-log scales. The de¬ 
pendencies obtained at different temperatures in a range 
from 6.25-14.00 Kelvin at zero magnetic field. At high 
applied currents the dependencies are in agreement with 
the power law V ~ where coefficient a depends on 
temperature. 

Fitting these curves to the power-law I-V function 
yields the temperature dependence of a(T). Figure 0 - 
presents the results of the fitting. The coefficient a varies 































5 


from about 7 at 6.25K to 1 at at 12K indicating a BKT 
transition at T c «7.5 K corresponding to a = 3. We 
note that in accordandce with the thermal conductance 
per unit area of the sample g e - p h ~ 1 • (T/20) 4 W/Kcm 2 
m the electron overheating at temperatures T=7-10K 
near the BKT transition is below 0.1 K at the highest 
currents applied to the sample. Thus contributions of 
thermal effects to the observed nonlinearity are negligi¬ 
bly small. 


NUMERICAL SIMULATIONS AND 
BOLOMETRIC CALIBRATION. 

The goal of the simulation is to use dc I-V character¬ 
istic of the sample to reproduce the I-V characteristics 
under the MW excitation. Figure [7] presents the approx¬ 
imation of the microwave circuit used in experiments. 
The circuit contains a 50 Ohms terminal resistance Rt 
which provides a broadband coupling of the circuit with 
microwaves. The nonlinear resistance R(t) describes the 
sample response. The contact resistance R c couples mi¬ 
crowaves with the sample.To isolate the dc and the mi¬ 
crowave circuits a 47 pF capacitor is added in series with 
the sample. The AC resistance of the capacitor is in¬ 
cluded in the R c . 

The simulation is based on the assumption that the 
nonlinear response to EM radiation follows the dc I-V 
characteristics. The response of the circuit is described 
by the following set of equations: 


V(t) = IsRc + Vs 

(3a) 

I(t) = I 0 + I s (t) 

(3b) 

Vo + V s (t) = V DC (I(t)) 

(3c) 


Eq.([3^) describes voltage in the microwave (AC) cir¬ 
cuit. The AC voltage V(t) at the end of the coax is the 
sum of the voltage applied to the contact resistance IsR c 
and the AC voltage applied to the sample Vs- where Is is 
the AC current through the sample and the resistance R c . 
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FIG. 6: (a) Shown here are the dc I-V characteristics for Sam¬ 
ple A to which a V ~ I a fitting has been done at high cur¬ 
rents, where the curves display a clear power-law form even at 
low temperatures, (b) Analysis of the temperature-dependent 
exponent a(T) demonstrates a BKT transition at T ~ 7-8K. 


Eq. ([sJd) describe the total current I ( t ) through the sam¬ 
ple, which is sum of AC (Is) and dc (Iq) currents. The 
third equation relates the voltage applied to the sample 
to the total current, using the dc V-I characteristic ob¬ 
tained in the experiment. 

The equations yield AC current Is ( t ) and voltage Vs (t) 
at given microwave voltage V(t) = V u cos(u)t ), contact 
resistance R c and dc current Jo- The time average of the 
AC response yields contributions to dc current 6Iq =< 
Is(t) > and voltage <5Vo =< Vs(t) > and makes a V-I 
dependence in the presence of the microwave radiation. 
The V-I characteristic depends on the contact resistance 
R c . A comparison between simulated and the actual V-I 
characteristics yields the contact resistance R c at given 
MW voltage (power P sim ). The tables below give the 
fitting parameter R c at different powers P sim . 


Figure 1. 


Psim(^W) 

R c (n) 

79.4 

240 

31.6 

200 

12.6 

180 


Figure 2. 


P S im(dbm) 

Rcfr) 

-40 to -15 

200 


Figure 3. 


Psim(^W) 

Rc{n) 

Pexp(nW) 

1.88 

61 

1.88 

3.74 

61 

3.74 

7.47 

61 

7.47 

0.5 

65 

13.0 

3.55 

65 

54.0 


We note that for all Figures except Fig.3 P S i m = P exp , 
where P exp is a calibrated power at the end of the coaxial 
line. 

Fig. 3 shows a reduction of the upper limit of the dc 
bias for the simulations at high MW powers. This is re¬ 
lated to the limited range of currents in the experimental 
dc V-I dependences used for the numerical simulation: 



FIG. 7: (color online) Shown here is the microwave circuit, 
which is used in the simulation to obtain the voltage and cur¬ 
rent in the sample. The circuit contains a terminal resistor 
(Rt), contact resistance (R c , a fitting parameter), and the 
sample itself (R(t), yellow). A dc bias Io is applied to the 
sample through dedicated leads equipped with low pass fil- 
ters(not shown). In addition the dc and microwave circuits 
are separated by 47 pF capacitor, which has been included in 

Rc 
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FIG. 8: Ratio between power P exp supplied to the end of 
the coaxial line at low temperatures and the calibrated power 
Pgen at the output of microwave generator vs frequency. The 
straight line is the approximation used to normalize the non¬ 
linear response shown in Fig. 5 of the main text. 


Mman4oi]- In accordance with Eqs. §3, c)the simula¬ 
tion calculates the voltage Vo + Vs = Vdc{I o + Xs(i)) us¬ 
ing the dc V-I dependence Vdc(I)- To obtain the voltage 
V the value Io + Is must be less I ma x- In Fig.3 it reduces 
progressively the range of the dc biases Iq at which the 
simulation works at higher MW currents (powers). 


Bolometric calibration of MW power. 

The in-situ bolometric calibration of the microwave 
power applied to the end of the coax is based on the 
following measurement. In experiment the resistor Rt = 
50fl terminating the coaxial line (see Fig[7f is thermally 
anchored to a copper sample holder. A microwave power 
P exp , applied to the end of the coax, heats the resistor 
Rt ■ It decreases the amount of the dc electric power, 
Pft , supplied by a temperature control unit to dc heater, 
which is attached to the sample holder to maintain the 
sample temperature. The change of the dc heating power 
A Ph = A (V£ c /Rh) is found by measuring voltage Vd c 
applied to the heater (resistance I?h,=100 f2) with and 
without MW radiation. The bolometric calibration of 
the MW power P exp is based on the conservation of the 
total energy supplied to the stage: A Ph = P exp - 

FigjH] shows the ratio of the calibrated power P exp to 
the power P gen at the output of the MW generator at dif¬ 
ferent frequencies. The calibrated power P exp decreases 
by about 10 times at high frequencies. The power de¬ 
crease is found to be in agreement (within 40% devi¬ 
ations) with total microwave losses in coaxial lines be¬ 
tween the MW generator and the sample, which were 
measured separately in transmission and reflection ex¬ 
periments. This agreement verifies the expected inde¬ 
pendence of the power losses in the terminal resistance 
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FIG. 9: Schematic for microwave reflection measurements 
with signals in the range of 1-7 GHz. The broadband split¬ 
ter splits microwaves into two channels. The reference signal 
goes through a directional coupler to a MW detector attached 
to port [1]. The incident signal goes down through another 
directional coupler and the probe, attached at port [2], to 
a sample stage attached to port [3]. The incident signal is 
reflected back through the probe to the MW detector (port 
[1]), where the interference occurs between the reflected and 
reference signals. 


Rt on the frequency. 

The bolometric calibration is done in the presence of 
the superconducting samples. The fact that the decrease 
in supplied MW power is largely due to MW losses in 
transmission lines indicates also a relatively weak effect 
of the superconducting samples attached to the end of 
the coaxial line on the overall microwave reflection from 
the sample stage. It is in accord with measurements of 
the reflection coefficient from the sample stage presented 
in the next section. The bolometric calibration allows 
computing the microwave voltage at the end of the coax: 
VL = {Pexp -Rt) 1 / 2 - The voltage is used in the numerical 
simulations. 


MICROWAVE COUPLING TO THE SAMPLES: 

REFLECTION MEASUREMENTS. 

In the low frequency domain the nonlinear behavior 
(and the contact resistance R c ) is found to be nearly fre¬ 
quency independent taking into account the bolometric 
calibration of the power applied to the end of the coaxial 
line. At high frequencies the response does not follow the 
I — V curves and the extraction of the contact resistance 
from the comparison with the numerical simulations is 
not possible. To evaluate the coupling of the EM radia¬ 
tion to the sample at high frequencies broad band mea¬ 
surements of a reflected microwave power were performed 
at two temperatures T n and T s corresponding to normal 
and superconducting states of the samples. Since in the 
normal state the superconductivity is absent the differ¬ 
ence between two reflected signals indicate the strength 
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of the coupling of the MW radiation to the superconduct¬ 
ing condensate at different frequencies. The microwave 
broad band measurements (1 to 7 GHz) demonstrate that 
microwaves are quite uniformly coupled to the supercon¬ 
ducting condensate in the studied frequency range indi¬ 
cating that frequency variations of the contact resistance 
R c are not significant for the observed strong reduction of 
the nonlinear response at high drive frequencies. Details 
of these experiments are presented below. 

The setup of the reflection experiments is shown in 
Fig|9] The microwave power supplied by a microwave 
generator splits between two channels. The reference 
(upper) channel provides a reference signal V c supplied 
through a broad band directional coupler (-10dB) to a 
broad band microwave detector, which is attached to port 

[1] , The sample (lower) channel supplies microwave ra¬ 
diation to another directional coupler (-16dB), which di¬ 
rects the radiation to the sample attached to port [3] of 
a semi-rigid coaxial line (probe) connected to the port 

[2] of the coupler. Reflected from the sample signal V r is 
guided by the same coax back to the detector, where it in¬ 
terferes with the reference signal V c . Due to a significant 
difference in lengths of the reference and sample chan¬ 
nels variations of the microwave wavelength (frequency) 
yield oscillating interference signal on the detector. At 
a small microwave amplitude the detected signal V det is 
proportional to microwave power P, in other words, to 
the square of the microwave voltage: 


V det ~ (V c + V r ) 2 (4a) 

^V 2 + 2-V c -V r -cos (j) + V 2 (4b) 

~ v 2 + 2 ■ R c • 7 2 rv m ■ cos </> + V 2 (4c) 


where these V c and V r are the reference and reflected 
voltages, considered on the phasor diagram and (f> is the 
phase difference between V c and V r signals. In the last 
equation the reflected signal V r = 7 2 TV)„ was substituted 
with the incident voltage Vt n at port [ 2 ], using the reflec¬ 
tion coefficient T at the bottom of the probe (port [3]) 
and the microwave losses 7 in the coaxial line between 
ports [2] and [3]. 

Terminating the port [3] by zero resistance (termi¬ 
nal impedance Z=0) one can measure the interference 
pattern corresponding to the complete reflection of mi¬ 
crowaves from the port [3]. Figure 10 a) presents the 
detector signal V de t corresponding to the complete mi¬ 
crowave reflection (Z=0) as function of microwave fre¬ 
quency. Fast oscillations corresponds to the interference 
between V c and V r , while the slow variations of the signal 
corresponds mostly to variations of the reference signal 
V c . Contributions of the reflected signal V 2 to the mi¬ 
crowave background are less than 10 % and are neglected 
below. To separate the interference term (fast oscilla¬ 
tions) from the background a FFT filter was applied re¬ 
sulting in the thick curve presented in FigJToJa) . A sub- 




FIG. 10: The above graphs demonstrate the procedure for ex¬ 
tracting the amplitude of the interference term between signal 
reflected from port [3] and the reference signal V c . The fig¬ 
ure presents the signal reflected from port [3] terminated by 
zero impedance (Z = 0). The same procedure is used for 
the signal reflected by the sample, (a) The full signal (green) 
is smoothed using a FFT filter (black), (b) Subtracting the 
smoothed signal from the full signal leaves only the interfer¬ 
ence term with frequency dependent amplitude shown in (c). 


traction of the background from the detector response 
yields the interference term 2 • V c ■ 7 2 rV) rl • cos <fi shown 
in Figfl0|b). In the case of the complete reflection the 
reflection coefficient To = —1. The corresponding am¬ 
plitude of the interference term 27 2 |r 0 | T 4 Fi,i is shown in 
FigJlOjc). A decrease of the interference amplitude be¬ 
low 2 GHz is related to the decrease of V c and V) due to 
a properties of the couplers, which are designed to work 
in the range between 2 and 7 GHz. Decrease of the am¬ 
plitude at high frequencies is related to microwave losses 
in the coaxial cable (y 2 ). 

When a sample attached to the port [3] the amplitude 
of the interference term is 2y 2 T V c Vi n , where T is the coef¬ 
ficient describing the microwave reflection from the sam¬ 
ple. One can see that the ratio between this interference 
term and the amplitude of the interference term corre¬ 
sponding to Z=0 yields the normalized interference term 
r • cos((()), which is proportional to the magnitude of the 
reflection coefficient from the sample. The phase tf> con- 



















































































Reflection at T=4.6K (a) 
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FIG. 11: (a) The frequency dependence of the reflection coef¬ 
ficient for Sample C at T =4.6K is found by dividing the inter¬ 
ference signal from the sample by the amplitude of the inter¬ 
ference signal corresponding to the perfect reflection (Z = 0 
terminal resistance shown in Figure 4c). (b) Curve 1 presents 
the change in the reflection coefficient from the sample stage 
attached to port [3], when the sample undergoes transition 
from the normal (T =13.6K) to superconducting (T =4.6K) 
states. Curve 2 presents the change of the reflection from the 
stage without the sample in the same temperature interval. 


tains a contribution related to the phase of the reflection 
coefficient. Thus the product T • cos(4>) contains both 
the amplitude and the phase of the reflection coefficient. 
Below we will name the product as reflection coefficient. 
Fig- 0 a) presents the reflection coefficient T • cos(</>) of 
Sample C as a function of the frequency. Sample C has 
the same structure as Sample B. 

The reflection coefficient depends on the impedance Z 
of the microwave circuit containing the sample. A low fre¬ 
quency approximation of the microwave circuit is shown 
in FigJT] The impedance depends on contributions from 
the 50 Ohm terminal resistor Rt, contact resistance R c 
and the sample. At high frequencies the impedance may 
contains contributions related to a geometry of the sam¬ 
ple holder, in particular, to finite lengths of the MW feed¬ 
ing leads. The change of the microwave reflection due to 
the superconductivity depends on the coupling between 
the sample and microwaves. Stronger coupling makes 
larger temperature variations of the reflection coefficient 


near the superconducting transition. The difference be¬ 
tween reflection coefficients in the superconducting and 
normal states is proportional to the superconducting cur¬ 
rent I s in the sample. The difference contains contri¬ 
butions related to variations of both amplitude and the 
phase of the reflection coefficient. Fig. |TT|[b) presents the 
difference A(T-cos(</>)) between two reflection coefficients 
obtained in superconducting (T=4.6K) state and normal 
(T=13.6K) states as a function of the frequency.The data 
show relatively weak and quite uniform variation of the 
difference A(T • cos(</>)) with the frequency indicating a 
uniform microwave coupling to the superconducting con¬ 
densate at different frequencies. In addition Fig. &) 
shows the temperature variation of the reflection from the 
stage without the sample. The trace indicates that the 
temperature variations of the reflection in the microwave 
setup itself are negligibly small [251 . 

To evaluate more quantitatively the microwave cou¬ 
pling we note that the current through the system is 
related to the circuit impedance Z: I u = V u /Z, where 
Vu is the microwave voltage applied to the circuit (port 
[3]). At small impedance variations the current changes 
are 

A/ w = A(^) = -^-Ar (5) 

, where we have used the relation between the impedance 
Z and the reflection coefficient F: ZfZ^ = (l + r)/(l —T) 
and the relation between incident calibrated voltage "/Vi n 
at port [2] and the voltage V u applied to the sample: 
Vu = lVin{ 1+ F). Z o = 50n is impedance of the coaxial 
line. 

Near the superconducting transition the impedance 
varies mostly due to the superconducting contributions. 
To evaluate the strength of the superconducting current 
I s we substitute the temperature change of the current 
A /„ in Eq. © between the normal and superconducting 
states by the superconducting current I s . In the case of 
the ’’parallel” electrical connection of the superconduct¬ 
ing layer to the end of the coaxial line the substitution 
is exact. For the circuit presented in Fig. [7] the relation 
between two currents is I s = (1 + R c /Rn)AI u , where 
Rn ~120 Ohm is the normal resistance of the sample. 
The contact resistance R c can be evaluated from the tem¬ 
perature variation of the reflection coefficient F. Indeed 
if the resistance R c is large in a comparison with the ter¬ 
minal resistance Rt one should expect a small tempera¬ 
ture variation of the impedance Z, since in both normal 
and superconducting states the current through the re¬ 
sistance R c and the sample is small. Figfllja) shows that 
the reflected power P r = T 2 {v)Pi n is less than 50% in the 
studied frequency range. It indicates that the impedance 
of the sample stage Z is quite close to the coax impedance 
Zq: Z rs Zq. Neglecting the difference one can related 
small variations of the reflection coefficient with small 
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FIG. 12: (a)Dependence of the microwave induced variation 
of the resistance A R on the square of microwave current, 
which was estimated in accordance with Eq.|8]| at different 
frequencies as labeled; (b) Frequency dependence of the non¬ 
linear coefficient k(u)=AR/ obtained in the limit of small 
MW currents and dc biases. Solid line presents a continu¬ 
ous frequency sweep. Red open circles are obtained from the 
dependencies presented in Fig |12[ a) using Eq.(|8|. All depen¬ 
dencies are obtained at T=8K. The insert shows the temper¬ 
ature dependence of the sheet resistance per square of Sample 
C. 

variations of the impedance: 

AF= (Z Z l A Z% * AZ/2Z 0 (6) 

(z + z Q y 

For the circuit shown in Fig. [7] the temperature vari¬ 
ation of the impedance between the normal (Rn) and 
superconducting (R$ ~ 0) states is 

A Z = R n ^ Rt _- R J n (7) 

Rt + Rc Rt + Rc + Rn 

Eqs. ([6]) and ([?]) yield that at Rt = 50 Ohm and 
AG > 0.1, which corresponds to Fig|TT|b), the contact 
resistance R c is below 70 Ohms in the studied frequency 
range. Thus the relation I s = A I u underestimates the 
actual superconducting current I s = (1 + R c / Rn)AI u 
by about 60%. This estimation error is much smaller 
than the observed frequency variations of the nonlinear 
response and is neglected below. 

Experiments indicate that at small currents the non¬ 
linear response is proportional to the applied power P 
and, therefore, to the square of the superconducting cur¬ 
rent I s . Eq.([5]) leads to the following evaluation of the 
nonlinear response at a small current: 

A R = m 2 s « k • ( 7 Ar) 2 • ^ (8) 

Zo 

, where k(^) is a frequency dependent nonlinear coeffi¬ 
cient and Pi n = V. 2 n j Z 0 is the calibrated power of the 
incident radiation at port [ 2 ], 

Figl^a) presents the dependence of the microwave 
induced change of the resistance on the square of the 


superconducting current at different frequencies as la¬ 
beled. The response, which is similar to the one ob¬ 
tained on other samples, shows the strong decrease in 
the magnitude at high driving frequencies. Figfl2^b) 
presents the frequency dependence of the nonlinear coef¬ 
ficient n(v)=AR // 2 . The dependence is obtained in the 
regime of a small power Pi n at which the nonlinear re¬ 
sponse is proportional to the square of the superconduct¬ 
ing current I s . The figure demonstrates again the strong 
reduction of the response to the applied microwave cur¬ 
rent at high driving frequencies. 

At low frequency v ~1 GHz the obtained coefficient of 
nonlinearity k(v) is in accord with the one obtained in 
dc domain experiments in a similar sample (sample S2 
in [T5]). The sample S2 has a slightly higher ( ST C ~1 
K) transition temperature. In the dc experiments the 
nonlinear coefficient was obtained from the relation E = 
PoJdc + 7 nlJcLc between the electric field E and applied 
current density Jdc = I/W, where W = 200/r is the 
width of samples and "/nl is the nonlinear coefficient. 
The relation can be rewritten as V = EL = Rq ■ I + 
”/nlL(I/W) 3 yielding k = 3jnlL/W 3 , where L = 800/i 
is the length of samples. At T=9K (corresponding to 
T= 8 K for Sample C) the nonlinear coefficient 'Jnl ~5 
Vim 2 /A 2 (see Fig.3a in [IS]). It yields «S 2 = 1 - 5 T 0 9 f 2 /A 2 . 
This value is in a quantitative agreement with the one 
shown in Fig | 12[ 3 at low frequencies. The correspondence 
between the nonlinear coefficients obtained by different 
experimental methods demonstrates that the reflection 
method provide a reasonable evaluation of the applied 
microwave current. 

In conclusion, the presented reflection measurements 
indicate that the microwave radiation is quite uniformly 
coupled to the superconducting electrons at different fre¬ 
quencies. Analyzed in terms of the applied superconduct¬ 
ing current the nonlinear response demonstrates signifi¬ 
cant reduction at high driving frequencies. 


* Corresponding author: vitkalov@sci.ccny.cuny.edu 
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